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Abstract 

The present paper reports interesting new features that wormhole solutions 
in the scalar field gravity theory have. To demonstrate these, we obtain, by us- 
ing a slightly modified form of the Matos-Nunez algorithm, an extended class 
of asymptotically flat wormhole solutions belonging to Einstein minimally cou- 
pled scalar field theory. Generally, solutions in these theories do not represent 
traversable wormholes due to the occurrence of curvature singularities. How- 
ever, the Ellis I solution of the Einstein minimally coupled theory, when Wick 
rotated, yields Ellis class III solution, the latter representing a singularity-free 
traversable wormhole. We see that Ellis I and III are not essentially independent 
solutions. The Wick rotated seed solutions, extended by the algorithm, contain 
two new parameters a and 5. The effect of the parameter a on the geodesic mo- 
tion of test particles reveals some remarkable features. By arguing for Sagnac 
effect in the extended Wick rotated solution, we find that the parameter a can 
indeed be interpreted as a rotation parameter of the wormhole. The analyses 
reported here have wider applicability in that they can very well be adopted in 
other theories, including in the string theory. 



I. Introduction 

Recently, there is a revival of interest in the scalar field gravity theories in- 
cluding the Brans-Dicke Theory due principally to the following reasons: Such 
theories occur naturally in the low energy limit of the effective string theory 
in four dimensions or the Kaluza-Klein theory. It is found to be consistent 
not only with the weak field solar system tests but also with the recent cos- 
mological observations. Moreover, the theory accommodates Mach's principle. 
(It is known that Einstein's General Relativity can not accommodate Mach's 
principle satisfactorily). All these information are well known. 
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A less well known yet an important arena where Brans-Dicke Theory has 
found immense applications is the field of wormhole physics, a field recently 
re-activated by the seminal work of Morris, Thorne and Yurtsever [1]. Concep- 
tual predecessors of modern day's wormholes could be traced to the geometry of 
Flamm paraboloid, Wheeler's concept of "charge without charge" , Klein bottle 
or the Einstein-Rosen bridge model of a charged particle [2]. Wormholes are 
topological handles that connect two distant regions of space. These objects 
are invoked in the investigations of problems ranging from local to cosmological 
scales, not to mention the possibility of using these objects as a means of inter- 
stellar travel [1] . Wormholes require for their construction what is called "exotic 
matter" - matter that violates some or all of the known energy conditions, the 
weakest being the averaged null energy condition. Such matters are known to 
arise in quantum effects (Casimir effect, for example). However, the strongest 
theoretical justification for the existence of exotic matter comes from the notion 
of dark energy or phantom energy that are necessary to explain the present ac- 
celeration of the Universe. Some classical fields can be conceived to play the role 
of exotic matter. They are known to occur, for instance, in the R + R 2 theories 
[3], Visser's thin shell geometries [4] and, of course, in scalar-tensor theories [5] 
of which Brans-Dicke theory is a prototype. There are several other situations 
where the energy conditions could be violated [6]. 

Brans-Dicke theory describes gravitation through a metric tensor ( g^ u ) and 
a massless scalar field ( <fi) . The Brans-Dicke action for the coupling parameter 
lo = — 1 can be obtained in the Jordan frame from the vacuum linear string 
theory in the low energy limit. The action can be conformally rescaled into 
what is known as the Einstein frame action in which the scalar field couples 
minimally to gravity. The last is referred to as the Einstein minimally coupled 
scalar field theory. Several static (mostly nontraversable) wormhole solutions 
in Einstein minimally coupled scalar field theory and Brans-Dicke theory have 
been widely investigated in the literature [7]. However, to our knowledge, exact 
rotating wormhole solutions are relatively scarce except a recent one in Einstein 
minimally coupled scalar field theory discussed by Matos and Nunez [8] . In this 
context, we recall the well known fact that the formal independent solutions 
of Brans-Dicke theory are not unique. (Of course, the black hole solution is 
unique for which the Brans-Dicke or minimal scalar field is trivial in virtue 
of the so called "no scalar hair" theorem.) Four classes of static Brans-Dicke 
theory solutions were derived by Brans [9] himself way back in 1962, and the 
corresponding four classes of Einstein minimally coupled field theory solutions 
are also known [10]. But recently it has been shown that only two of the four 
classes of Brans' solutions are independent [11]; the other two can be derived 
from them. However, although all the original four classes of Brans' or Einstein 
minimally coupled solutions are important in their own right, we shall here 
consider, for illustrative purposes, only one of them (Ellis I) as seed solution. 
The same procedure can be easily adopted in other three classes. 

The general motivation in the present paper is the following: To properly 
frame an algorithm for generating singularity-free asymptotically flat rotating 
wormhole solutions from the Ellis seed solutions and to investigate the role of 
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new parameters in the extended solutions. The analyses also answer a certain 
long standing query on the wormhole solutions in the Brans-Dicke theory. 

In this article, using a slightly modified algorithm of Matos and Nunez [8] , 
we shall provide a method for generating wormhole solutions from the known 
static seed solutions belonging to Einstein minimally coupled scalar field theory. 
The solutions can then be transferred to those of Brans-Dicke theory via inverse 
Dickc transformations. For illustration of the method, only Ellis I seed solution 
is considered here, others are left out because they can be dealt with similarly. 
The Brans-Dicke solutions can be further rephrased as solutions of the vacuum 
4-dimcnsional low energy string theory ( u> = — 1) and the section II shows how 
to do that. In sections II- V, we shall analyze and compare the behavior Ellis III 
and the Wick rotated Ellis I solution pointing out certain interesting differences 
in these geometries. In section VI, the study of the geodesic motion in the 
extended Wick rotated Ellis I solution reveals the meanings of the Matos-Niinez 
parameter a and section VII shows, via consideration of the Sagnac effect, that 
a can indeed be accepted as a rotation parameter. Finally in section VIII, 
we shall summarize the results. Throughout the article, we take the signature 
(—,+,+,+) and units such that 8irG = c = 1, unless restored specifically. Greek 
indices run from to 3 while Roman indices run from 1 to 3. 



Let us start from the 4-dimcnsional, low energy effective action of heterotic 
string theory compactified on a 6-torus. The tree level string action, keeping 
only linear terms in the string tension a' and in the curvature R, takes the 
following form in the matter free region ( S matter = 0): 



where g^ v is the string metric and <J> is the dilaton field. Note that the zero 
values of other matter fields do not impose any additional constraints either on 
the metric or on the dilation [12]. Under the substitution e~ 2 * = <j>, the above 
action reduces to the BD action 



in which the BD coupling parameter u> = —1. This particular value is actually 
model independent and it actually arises due to the target space duality. It 
should be noted that the BD action has a conformal invariance characterized 
by a constant gauge parameter £ [13]. Arbitrary values of £ can actually lead 
to a shift from the value u> = — 1, but we fix this ambiguity by choosing £ = 0. 
Under a further substitution 



II. The action, ansatz and the algorithm 
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in which we have introduced, on purpose, a constant parameter a that can have 
any sign. Then the action (2) goes into the form of EMS action 

S E ms = J d 4 x^ [R + • (4) 

The EMS field equations are given by 

R^ v = -ol<p^ v (5) 

K = ( 6 ) 

We shall choose a = +1, ip = <p{l) in what follows. The negative sign on the 
right hand side of Eq.(5) implies that the source stresses violate some energy 
conditions. The ansatz we take is the following: 

ds 2 = - f (I) (dt + a cos Odipf + f' 1 (I) [dl 2 + {l 2 + l 2 ) (dO 2 + sin 2 9d^ 2 )] , (7) 

where lo is an arbitrary constant, the constant a has been interpreted in [8] as a 
rotational parameter of the wormhole. We call it the Matos-Nunez parameter. 
The ansatz in (7) is actually a subclass of the more general class of stationary 
metrics given by [14,15]: 

ds 2 = -/ (dt - uj l dx i ) 2 + f- 1 h i:j dx i dx j (8) 

where the metric function /, the vector potential and the reduced metric hij 
depend only on space coordinates x l . We shall see below that the parameter 
a can be so adjusted as to make a symmetric, traversable wormhole out of an 
asymmetric, nontraversable one. Note also that the replacement of a — > —a 
does not alter the field equations. 

The function f(l) of the ansatz (7) is then a solution of the field equations 
(5) and (6) if it satisfies the following: 



(I 2 + It) J 



a 2 f 2 



\f) (P + 'S) 

where the prime denotes differentiation with respect to I. 
Algorithm: 

Let /o = fo(l',p, q;a = 0) and ip = ip (l;p, q; a — 0) be a known seed solution 
set of the static configuration in which p, q are arbitrary constants in the solution 
interpreted as the mass and scalar charge of the configuration. Then the new 
generated (or extended) solution set (/, ip) is 

f(l;p,q;a) = — - T — ^, ip(l;p, q; a) = <p (11) 

a + nd f 
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where n is a natural number and the parcmctcrs p, q are specific to a given 
seed solution set (f ,ip ) while S is a free parameter allowed by the generated 
solution in the sense that it cancels out of the nonlinear field equations. The 
scalar field <p is remarkably given by the same static solution of the massless 
Klein-Gordon equation tp : £ = 0. The seed solution ( a = 0) following from 
Eqs.(9) and (10) gives S — 2pq. For the generated solution (a ^ 0), the value of 
d may be fixed either by the condition of asymptotic flatness or via the matching 
conditions at specified boundaries. Eq.(ll) is the algorithm we propose. This 
is similar to, but not quite the same as, the Matos-Nunez [8] algorithm. The 
difference is that they defined the free parameter as 5 = \f~D . The difficulty in 
this case is that, for our seed solution set (/ , f ) below, the field Eqs.(10) and 
(11) identically fix S 2 = D = giving / = which is obviously meaningless. The 
other difference is that we have introduced a real number n that now designates 
each seed solution / and likewise the corresponding new solution /. With the 
known parameters n, p and q plugged into the right side of Eq.(ll), the new 
solution set (f,ip) identically satisfies Eqs. (9) and (10). One also sees that the 
algorithm can be applied with the set (/, ip) as the new seed solution and the 
process can be indefinitely iterated to generate any number of new solutions. 
This is a notable generality of the algorithm. 

III. Ellis I solution and its geometry 

The study of the solutions of the Einstein minimally coupled scalar field 
system has a long history. Static spherically symmetric solutions have been 
independently discovered in different forms by many authors and their properties 
are well known [16,17]. We start from the following form of Class I solution, 
due to Buchdahl [18], of Einstein minimally coupled theory : 



ds 2 



-i rn 

1 2r 



1 + 



2r 



26 



dt 1 



2(l+/3) 



[dr 2 + r 2 d0 2 + r 2 sin 2 Bd^ 2 \ 



(12) 



Q 



2r 



1 + — 

1 T 2r 



(13) 



where m and /3 are two arbitrary constants. The same solution, in harmonic 
coordinates, has been obtained and analyzed also by Bronnikov [19]. 
The metric (13) can be expanded to give 



ds 2 



1 



2M 
r 



2M 2 



oil 



dt 2 



[dr 2 +r 2 d8 2 + r 2 sin 2 0# 2 ], 
from which one can read off the Keplerian mass 

M = m/3 



1 



2M 
r 



(14) 



(15) 
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The solution (12) describes all weak field solar system tests because it exactly 
reproduces the PPN parameters. For [3 = 1, it reduces to the Schwarzschild 
black hole solution in isotropic coordinates. For a = +1 and (i > 1, it represents 
a naked singularity. The metric is invariant in form under inversion of the radial 
coordinate r — ► ^- and we have two asymptotically flat regions (at r — and 
r = oo), the minimum area radius (throat) occurring at ro = ^ j3 + \/ (3 2 — 1 . 

Thus, real throat is guaranteed by the condition (3 2 > 1 which we might call here 
the wormhole condition. However, despite these facts, because of the occurrence 
of naked singularity at r = m/2, it is not traversable and so Visser [4] called it 



a "diseased" wormhole. For the choice a = +1, the quantity \J2{0 2 — 1) is real 
such that there is a real scalar charge a from Eq.(14) given by 



a 2m f3 2 - 1 

^ = 7 = -vVV (16) 

But, in this case, we have violated almost all energy conditions importing by 
hand a negative sign before the kinetic term in Eq.(5). Alternatively, we could 
have chosen a = — 1 in Eq.(13), which would then give an imaginary charge 
a. In both cases, however, we end up with the same equation = —tp ip v . 
There is absolutely no problem in accommodating an imaginary scalar charge 
in any configuration violating energy conditions [17,18]. 
The Ricci scalar R for the solution (12) is 

R _ 2mV(l-/3 2 ) 

(r-ro/2) 2 ( 2 -«(r + m/2) 2 ( 2 +« 1 ' 

which diverges at r = m/2 showing a curvature singularity there. For (3 > 2, the 
divergence in the Ricci scalar is removed, but then the metric becomes singular. 
However, metric singularity is often removable when one redefinines it in better 
coordinates and parameters. 

Using the coordinate transformation I = r + ^r, the solution (12) and (13) 
can be expressed as 

ds 2 = -f (l)dt 2 + -4tt [dl 2 + (I 2 - m 2 ) (d9 2 + sin 2 6dip 2 )} , 

/ow = ( i T^y, as) 




In this form, it is exactly the Ellis I solution that has been discussed also by 
Bronnikov and Shikin [20]. Eqs.(18), (19) identically satisfy the field Eqs.(9) 
and (10) for a = 0. This is our seed solution set but we still need to suitably 
redefine it because of the apprearance of the naked singularity at / = m The 
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throat lo appears at Z = r o + ^ = m P > m corresponding to r = r . Thus 
the minimum surface area now has a value 4nm 2 f3 2 . For this solution, p = m, 
q = /?, and it turns out that the seed solutions (18), (19) correspond to n = 4 
so that the generated EMS solution that identically satisfies the field equations 
(9) and (10) for a ^ are: 

/(Z; to, /?, a) = — — 2^2 ^(0 = y ln 

To achieve asymptotic flatness at both sides, that is, /(Z) — > 1 as Z — > ±oo, we 
note that /o(Z) — > 1 as Z — > ±oo . Therefore, we must fix 

c 2M ± V4M 2 - a 2 
S = 2 

In the above, we should retain only the positive sign before the square root. The 
reason is that, for a = 0, the negative sign gives 8 = implying / = which is 
meaningless. On the other hand, the positive root gives 8 = 2M and / = /o, 
as desired. Note that = 1 does not lead to Kerr black hole solution from the 
generated metric. Therefore, the latter does not represent a rotating black hole 
but might represent the spacetime of a rotating wormhole [8] . 

Let us now examine wormhole geometries in the static and generated solu- 
tions in the EMS. 

(a) Static seed case (a = 0): 

The first observation is that the metric functions in Eq.(18) diverge at the 
singularity Z = ±m as does the Ricci scalar. The next observation relates to 
the behavior of the area radius. It exhibits certain peculiar properties for the 
metric in (18) for (3 > 1. For the segments Z > m and Z < —to, we have the 

area radius pl(l) = \J f^il 2 - to 2 ). Then, the area Airp^{l) decreases from +oo 
at one asymptotic flat end to a minimum value p min = /°o('o) & t the throat 
Z = Zo = to/3, and then becomes asymptotically large, but not flat, at a radial 
point I — m. In the remaining segment, we have |Z| < to, and the area now has 

to be redefined as Pq 7 (Z) = \J Iq 1 \m 2 — l 2 \. It then decreases to zero at Z = —to 
(another throat, zero radius!) and then opens asymptotically out to — oo at the 
other asymptotic flat end. Though in the r— coordinate version, the metric is 
inversion symmetric, there are now two asymptotically flat, isometric universes 
with their own throats and are actually disconnected by the naked singularity 
at Z = to. They are also asymmetric around Z = due to the fact that the 
throat radii in the two universes are different. 

(b) Generated case (a ^ 0): 

Here again, the metric functions given by Eq.(20) diverge at Z = ±to. The 
throat of the rotating wormhole can be found from the roots of the equation for 
Z : 

a 2 [l{f 2 - 1) + mP((f 2 + 1)] + 4to/3(Z - to/3)(2to/3 + ^Am 2 (3 2 -a 2 ) = (22) 
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TO 



Z + TO 



(20) 



(21) 



They can be computed only numerically for given values of the parameters to 
and (3. However, the area radius p(l) = ^// _1 (Z 2 — m 2 ) for the solution / shows 
that, for a / 0, the area jumps to infinity for (3 > 1 at I = ±m but flares 
out asymptotically to ±oo on both sides. One now has three disconnected uni- 
verses, that is, a one-sided asymptotically flat universe, a both-sided non-flat 
but asymptotically large "sandwich universe" and another one-sided asymptot- 
ically flat universe. For the extreme case a = 2m(3, the picture is the same 
but the behavior of p(l) is symmetric around I = 0. Thus, the wormhole is not 
traversable, be it seed solution (18) or extended solution (20). A natural ques- 
tion arises if they be made traversable by removing the singularity manifested 
in the infinite jump in the area at I = ±m and in the curvature. We shall 
consider this question now. 

IV. Ellis III solution via Wick rotation 

One procedure to remove the aforementioned singularities is to analytically 
continue the Ellis I solutions (fo:fo) by means of Wick rotation of the param- 
eters while maintaining the real numerical value of the throat radius. In the 
solution set (/ ,</? ), we choose 

to -> -im, [3 -> i(3 (23) 

so that l = m (3 is invariant in sign and magnitude. 
(a) a = 0: 

Then the metric resulting from the seed Eq.(18) is and it is our redefined 
seed solution: 

ds 2 = -f' {l)dt 2 + -J— [dl 2 + (l 2 + m 2 ) (d9 2 + sin 2 dd^ 2 )} (24) 



io(0 = cx p 



—2(3 arccot 



(25) 



Po(0 = 



1 + /3 2 



arccot ( — I (26) 



/// 



This is no new solution but is just the Ellis III solution [19] which can be 
obtained in the original form by using the relation [20] 

arccot(a;) + arctan(x) = + — ;x>0 (27) 

= -|;*<o (28) 

and the function on the left shows a finite jump (of magnitude n) at x = 0. 
Thus, we get from Eqs.(25) and (26) two branches, the +ve sign corresponds to 
the side I > and the —ve sign to I < [22]: 

f E l lis (l) = cxp[-2/?{±5 - arctan(-)}] (29) 
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Ellis fj\ 



V2 yj 1 + /3 2 



± — — arctan( — ) 



(30) 



We might study the solutions (30) and (31) per se, or equivalently, study the 
two restricted branches taken together, while allowing for a discontinuity at the 
origin I = 0. Otherwise, we might disregard (30), (31) and treat each of the 
± set in Eqs.(34), (35) as independently derived exact solution valid in the un- 
restricted range of Z with no discontinuity at I = 0. The two alternatives do 
not appear quite the same. In fact, each of the individual branch represents 
a geodesically complete, asymptotically flat traversable wormhole (termed as 
"drainholcs" by Ellis) having different masses, one positive and the other nega- 
tive, on two sides respectively. The known Ellis III solution is the +ve branch 
which is continuous over the entire interval Z G (— oo, +oo). The — ve branch is 
also equally good. What we have shown here is that the Ellis solutions I and 
III are not independent solutions of the Einstein minimally coupled scalar field 
theory as one can be obtained from the other. 

It is of interest to compare the behaviors of the Ellis III solutions (34) 
with the Wick rotated Ellis I solutions (30): (i) The Ellis III metric function 
fo+ is (l) -> 1 as Z -> +oo but ,f^l Us (l) -» e- 27r/3 as I -> -oo. These two limits 
correspond to a Schwarzschild mass M at one mouth and -Me"' i at the other. 
There is no discontinuity at the origin because fof hs (l) — > e -71 "^ as I — > ±0. 
In the solution (30), on the other hand, there is a discontinuity at the origin 
because /o(Z) — > e ±7r ' 3 as Z — > ±0 while there is no asymptotic mass jump since 
fo(l) — > 1 as I — > ±oo. The curvature scalars for both (30) and (34) are formally 
the same and given by 



R'n 



jjEllis 
K 0+ 



2m 2 (l±/3 2 ) 
(Z 2 +to 2 ) 2 

2to 2 (1 + /? 2 ) 
(I 2 + m 



exp 



■ exp 



-2/3 arccot 



-2/3{— — arctan 



(31) 
(32) 



which go to zero as Z — ► ±oo. That means, the spacetime is flat on two sides for 
both the solutions. Next, we verify what happens to these scalars at the singular 
coordinate radius (r = to/2) that has now been shifted to the origin I = r— ^ = 
0. (ii) The Ellis curvature scalar R^ ls -» - %±±fl e -*P as I -» ±0 whereas the 
curvature scalar R' exhibits a finite jump from — - e^* 13 to — - e+7T/3 
as I -» ±0. (iii) The area radius p^ is (l) = y/ fo+ iEUis) (I 2 + m 2 ) -» toVP^ as 

Z — ► ±0 whereas Po(Z) = \J~f^~H} 2 + ™ 2 ) shows a finite jump from mV to 

mVe~^P as Z — > ±0. These show that while Ellis wormhole (34) is traversable, 
but jumps at the origin in the Wick rotated solution (30) prevent traversability. 
The behaviors of (34) and (30) are different at the origin except that, for both 
the solutions, the throat appears at the same radius Zo = M = mf3. Similar 
considerations apply for the —ve branch. 

Ellis III wormholes (34) can be straightforwardly extended to the rotating 
form via the algotithm (11) and they would be likewise traversable, as has been 
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shown by Matos and Nunez [8]. Thus we do not deal with it here. Henceforth, 
we would rather concentrate on the Wick rotated Ellis I solution and ask: Can 
we somehow remove the discontinuities in (30)? That is exactly where the new 
paramter a comes in. Let us see what happens when a^O. 
(b) a £ : 

The minimum area radius of the extended solution is obtained from the 
equation -jjj- = where p'(l) — \//' -1 (Z 2 + m 2 ) is the area radius. From nu- 
merical study of the resulting equation, we find that the minimum area occurs 
at I < m(3 and it decreases with the increase of a for fixed values of m and (3. 
We also notice that the finite jump persists in the area radius of the extended 
solution f'(l; m, 13, a), 8 being still given by Eq.(21). Surprisingly however, when 
a = 2m(3, the area function p'(l) decreases from +oo to the minimum value at 
the throat, then increases to a finite value at I = 0, undergoes no jump at I = 
but passes continuously, though not with C 2 smoothness, across I = 0, on to 
— oo. The Ricci scalar R' for f'(l; m, (3, a) is given by 



, 8/3 (1 + (3 2 ) (2m(3 + x/4m 2 (3 2 ~ a 2 ) exp[2/?arccot(^)] 
R = 



m 



(2m (3 + \/ Am 2 (3 2 — a 2 ^j cxp[4/3 arccot(- 



(33) 



and it approaches the value ^sT^gfe?) as ' ~~ * that is, no jump in it. 

The area behavior shows that, for the extreme case a = 2m(3, we do have a 
traversable wormhole with a single metric covering both the asymptotically flat 
universes (I — > ±oo) connected by a finite wedge-like protrusion in the shape 
function at I = 0. This wedge prevents C 2 continuity across I = in the area 
function but sews up two exactly symmetrical asymptotically flat universes on 
both sides. The numerical values of the free parameters m and (3 can always 
be suitably controlled to make the tidal force humanly tolerable and the travel 
safer. 

For a ^ 2m(3, such a single coordinate chart is not possible as the area has 
a jump at I = 0. However, we can artificially circumvent this discontinuity 
and connect the two disjoint universes by multiple metric choices on different 
segments. We can get a cue for this construction from the static case. Consider 
the metric form (18) on one segment (AB) and the Wick rotated metric (24) on 
the other side (BC) so that the areas match at a radial point I = l\. The radius 
I = l\ is a root of the equation (area from right (AB) = area from left (BC)) 

(m 2 - I 2 ) x (j^^j ^ = (m 2 + I 2 ) x £zp[2/?arccot(^)] (34) 

By numerical computation, it turns out that < l\ < m such that the two 
otherwise disjoint universes, one represented by the branch AB and the other 
by BC, can be connected at B(l = l\). At the joining point B, there is continuity 
in the area function (again not C 2 ) and the tidal forces can be shown to be finite 
throughout the generator curve ABC. Exactly similar arguments hold in the 
rotating case. Branche AB belongs to the Wick rotated solution (/') while the 



10 



sector BC belongs to the original solution (/). Numerical calculations show 
that the matching occurs at either of the two points B{1\) or B(l 2 ) such that 
-m < h,h < m. 

Traversable wormholes can also be constructed by employing the "cut-and- 
paste" procedure [4]. One takes two copies of the static wormholes and joins 
them at a radius I — l b > l . The interface between the two copies will then 
be described by a thin shell of exotic matter. The shape functions on both 
sides will be symmetric. However, when rotation is introduced, numerical cal- 
culation shows that the throat radius decreases from the static value while the 
flaring out occurs faster. It is of some interest to note that Crisostomo and Olca 
[20,21] developed a Hamiltonian formalism to obtain the dynamics of a massive 
rotating thin shell in (2+1) dimensions. There, the matching conditions are un- 
derstood as continuity of the Hamiltonian functions for an ADM foliation of the 
metric. Of course, this procedure can be trivially extended to deal with axially 
symmetric solutions in (3+1) dimensions. For soliton solutions, see Ref.[22]. 

V. Extended Brans-Dicke I solution 

To obtain the rotating Brans-Dicke solution, we pursue the following steps: 
Note from Eq.(3) that 



2uo 



■ In , 



<f — In 



m -, V2(/? 2 -l) 
1r 



1 + — 

1 + 2r 



1 - 



2r 



1 + — 



V4(/3 2 -l)/(2^+3) 



(35) 



Now using the constraint from the Brans-Dicke field equations [9], viz., 

4(/3 2 -l) = -(2^ + 3)^, (36) 

where C, A are two new arbitrary constants and ui is the coupling parameter, 
we get 

c c 

(37) 



"1 - 


m " 

1r 


c 

X 


I — m 


C 
2A 


.1 + 


m 
2r - 




I + m 





The Eq.(34) can be rephrased in the familiar form [9] 

n2 W(7 



A 2 = (C+1) 2 -C(1 



)• 



(38) 



However, the minimum area condition 1 > 1 requires that the right hand side 
of Eq.(34) be positive. This is possible if either u < — | or A be imaginary. Let 
us first consider u) < — | so that the exponents are real. Then, the final step 
consists in using the relation = (j)^ 1 g tn y together with replacing (3 in the 
exponents in the g^ u by [7] 



(39) 
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This means we have the Brans-Dicke rotating class I solution for a; < — § 
follows: 



MS 



ds 2 



-h{l)dt 2 + / 2 (0 [dl 2 + (I 2 - to 2 ) (d9 2 + sin 2 6d^ 2 ) 
Ji (I; to, C, A, a) = f(l; m, (3, a)^ -1 

8mS[+ x (C + 2)] f * lHrJ 



l-\-m 



a 2 + AS 2 



l+m 



/ + TO 



(40) 



(41) 



h{i) 



f 2 {l\m,C,\a) = f 1 (l;m,(3,a) 



A5 2 



l—m 



l+m 



l — m 



l+m 



/ + TO 



C 
" 2A 



m = 



TO 



/ + TO 



(42) 



(43) 



It can be verified that the Brans-Dicke field equations again yield the expression 
(36). Using the relation I = r + 2p, it can be easily expressed in the familiar 
(i, r, 8, ip) coordinates with the value of 8 given by Eq.(21) in which /3 should have 
the value as in Eq.(36). For instance, when a = 0, we have S — ^(C + 2) and 
identifying ^ = B, one retrieves the static Brans-Dicke metric in the original 



notation: 



2(A-C-1) 



ds 2 = 



1-^ 



dt z 



B 

1 + — 

r 



1- f 

T+t 



[dr 



2 d6> 2 + r 2 sin 2 6»d^ 2 



1-S 

r_ 

i + f 



(44) 
(45) 



The solution (42) represents a naked singularity at r — B and the condition for 
the existence of a minimum area is [6,7] 



(C + l) 2 > A 2 . 



(46) 



For (3 > 1, and a = +1, the negative kinetic term in the field equations (5) 
shows that the energy density is negative violating at least the Weak Energy 
Condition. The solution (38) still does not represent a traversable rotating 
wormholc in the Einstein minimally coupled theory. One has to first make a 
change 

to — > im, A — > —iX (47) 

in (38) and (41) to get the Wick rotated counterpart. The next step is to use 
the relations (32), (33) obtaining two branches as obtained in Sec. IV. Either 
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of the branches would then represent rotating traversable wormholes in the 
Brans-Dicke theory for the range of coupling values u) < — §. Same classes of 
solutions will be obtained by alternative calculations with f3 2 > 1 and a = — 1 
(imaginary scalar charge). The above steps represent the basic scheme to be 
followed in other classes of solutions in the Einsten minimally coupled or Brans- 
Dicke theory, which we do not do here. 

The discussion about traversability of the wormholes in the Einstein mini- 
mally coupled theory can be transferred almost in verbatim to the Brans-Dicke 
theory once we use the crucial relation (37) connecting the parameters in both 
the theories. 

VI. Geodesic motion in the extended solution 

We shall consider, as an illustration, the class of solutions (24) generated 
from the Wick rotated seed solutions (18) of Einstein minimally coupled theory. 
The resulting solution is (now dropping primes) 



8m/?£exp 2/?arccot(-M 

f(l;m,(3,a) = ' , PL ' 48 

a 2 + 4(5 2 cxp[4/3arccot(i) 1 



p(0 



VV+/3 2 ) 



arccot I — ) (49) 



in 



so that the metric components g v \ are 



5oo = -/, (50) 

3ii = r\ (51) 

322 = f-\l 2 + ll), (52) 

333 = f - 1 (I 2 + l 2 ) sin 2 6 -fa 2 cos 2 6, (53) 

303 = 330 = -fa cos 6. (54) 

The four velocity is defined by 

dr^ 

U» = —,x» = (t,l,0,4>),dp = mods (55) 
dp 

in which p is the new affine parameter and mo is the invariant rest mass of the 
test particle. The geodesic equations are given by 



dUf,, 1 dg v \ x 



dp 2 dx» 



U V U X = 0, g vX U v U x = m 2 = e (56) 



Since the metric functions g v \ do not contain t and <j>, the corresponding mo- 
menta are conserved, that is, the fi — and p — 3 equations give, respectively 

U = -f.U° - faU 3 cos9 = k (57) 
U 3 = -. f aU° cos 9 + U 3 [.r 1 (l 2 + ll) sin 2 9 -fa 2 cos 2 6} = h (58) 
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where k and h are arbitrary constants. From the above, it follows that the 
"angular momentum" of the particle is 

^ 3 = fe-fcacosfl 
sin 9 

r 2 (l) ee Z^^ + lg). (60) 
The fj, = 2 component of the geodesic equation gives 

T" (W)!?) + i/«^°^ 3 sin^- (tf 3 )V(0 + /a 2 ] cos sin = 0. (61) 
dp \ dp j £ 

Instead of the fi = 1 equation, we take the second of Eq.(53) which gives 

e = -f[U° + U 3 acos9} 2 + /^(U 1 ) 2 + r 2 (l)[(U 2 ) 2 + (U 3 ) 2 sin 2 9]. (62) 
Looking at Eq.(58), we see that two solutions are possible. One is 

U 3 = => 9 = arccos(-^-) = const. => U 2 = (63) 
ka 

which means 9 can assume any constant value depending on the independent 
values of h, k and a. But such motions will only be radial since U 2 = and 
U 3 = 0. In other words the gravitating source acts like a radial sink\ The 
geodesic equation of motion (59) reduces to: 

f^j =ef(l;m,[3,a) + k 2 (64) 
This can be rewritten very succinctly as 

d^=2di- (65) 

The other solution of Eq.(58) is 9 = 0, which implies that the test particle 
motion is restricted to polar planes. However, we can always choose the pole 
perpendicular to this plane so that the angle tp varies in that plane with the 
particle motion. Moreover, from Eq.(56), we must have h = ka so that r 2 (l)U 3 = 
£ 7^ O.Then, we end up again with the same metric function but without the 
explicit appearance of a as an arbitrarily regulated free parameter. Thus, in the 
equation of motion, we have / = f(l; m, [3, h/k) so that a, which is a parameter 
of the gravitating source, is obtained from the motional characteristics like h 
and k of the test particle itself. This intriguing feature is somewhat analogous to 
the fact that the mass of the gravitating Sun can be determined from the motion 
of a test particle (planet) around it. As a special case, the parameter a can be 
so chosen as to completely mask the angular momentum £ of the test particle, 
that is, as h — > ka, there is a possibility that £ — > too. Physically, it is like 
choosing the parameter a to coincide with the orbital tp— angular momentum of 
the test particle. The equation of motion is again exactly the same as Eq.(59) 
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since U 2 — even though U 3 (^ 0) does not appear explicitly. This is due to 
the fact that, in the last term of Eq.(58), (U 3 ) 2 sin 2 6 — but the signature of 
orbiting (non-radial) test particle is manifest in the presence of h: 

— J = ef(l;m,(3,h/k) + k 2 (66) 



dp 

The turning points of the orbit will occur when ef = —k 2 and | ^ 0. From 
these conditions, we have the turning points occurring at 

l = l Q =mcot(^^J (67) 

where 

-2m(3e ± v / 4m 2 2 e 2 - k 4 a 2 

x ° = 2m • (68) 

Circular orbits will occur if ef = —k 2 and % = and they will be stable if 
^ < and unstable if ^ > 0. 

VII. Sagnac effect 

The best way to assess the effect of a non-zero U 3 (or dip ^ 0) and of the 
Matos-Nunez parameter a is through the Sagnac effect [23] analyzing the orbit 
in the plane 9 — 0. The effect stems from the basic physical fact that the 
round trip time of light around a closed contour, when the source is fixed on 
a turntable, depends on the angular velocity, say fi, of the turntable. Using 
special theory of relativity and assuming Qr -C c, one obtains the proper time 
8t s , when the two beams meet again at the starting point as 

St s £* -3-S (69) 

where S (= irr 2 ) is the projected area of the contour perpendicular to the axis 
of rotation. 

Without any loss of rigor, we take a —> —a for notational convenience al- 
though it is not mandatory. Suppose that the source/receiver of two oppositely 
directed light beams is moving along a circumference I = R = constant. Suit- 
ably placed mirrors reflect back to their origin both beams after a circular trip 
about the central rotating wormholc. (The motion is thus not geodesic or force 
free!). Let us further assume that the source/receiver is moving with uniform 
orbital angular speed ui n with respect to distant stars such that the rotation 
angle is 

<fio = ^ot- (70) 
Under these conditions, the metric becomes 

dr 2 =-/(i?;m,/3,a)[l-aw ] 2 ^ 2 (71) 
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The trajectory of a light ray is dr 2 = which gives 



[1 - au] 2 = 



(72) 



where u is the angular speed and where we have suspended the condition h = ka 
for photon orbits. The roots of the above equation coincide and is 



u>i± 



1 

a 



Therefore, the rotation angle for the light rays is 



If = U>l±t = U>1±- 



U) 



(73) 



(74) 



The first intersection of the world lines of the two light rays with the world line 
of the orbiting observer after emission at time t = occurs when 



so that 



ip + = ip + 2ir, ip 



U! 



(75) 



(76) 



where ± refer to co-rotating and counter-rotating beams respectively. Solving 
for ip , we get 

¥>o± = ± • (77) 

wi± - 

The proper time as measured by the orbiting observer is found from Eq. (66) by 
using dt = dip /oJo and integrating between ip 0+ and ip Q _. The final result is 
the Sagnac delay given by 



\6t„ 



= ^f{R;m,f3, a) 
= V7(RyrnJ3, a) 



1 — clluq 

1 — aujQ 
uj 



(^0+ - Vo-) 



■ w 



= y / .f(R;m,P,a)(4TTa). 



(78) 



This shows remarkably that the Sagnac delay depends only on the Matos- Nunez 
parameter a. Interestingly, the result is independent of w meaning that it is 
independent of the motional state of the source/receiver, be it static with respect 
to the frame of distant stars or moving with regard to it. When a = 0, there is 
no delay because the wormhole spacetime is then nonrotating (no turntable!). 
The above result supports the conclusion of Ref. [8] from an altogether different 
viewpoint that a can indeed be interpreted as a rotational parameter of the 
wormhole. 



VIII. Summary 
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Asymptotically flat rotating solutions are rather rare in the literature, be 
they wormholcs or naked singularities. The algorithm (11), together with some 
operations, provides a method to generate new traversable wormhole solutions 
in the Einstein minimally coupled and then in the Brans-Dicke theory. The 
present study opens up possibilities to explore in more detail new solutions in 
other theories too. For instance, the string solutions are just the Brans-Dicke 
solutions with u> = — 1. As we saw, the asymptotically flat wormhole solutions 
admit two arbitrary parameters a and 6. The Matos-Nuncz parameter a has 
been interpreted by the authors (Ref [8]) as a rotation parameter. Here we have 
shown how the parameters can be adjusted for obtaining traversability. 

The static Ellis I wormholes do not appear traversable due to the singularity 
as manifested in the behavior of the area function and curvature. To tackle 
this problem, we analytically continued it via Wick rotation and rederived El- 
lis III singularity free, asymptotically flat, traversable wormholes as one of its 
branches. This showed that the two classes of solutions are not independent, 
contrary to the general belief. The comparative features of the Wick rotated 
and Ellis III solutions are pointed out. In the extended Wick rotated solution, 
numerical graphics showed that the wormhole can be covered by a single metric 
with C° continuity. That is, the jumps can be sewed up at the origin for the 
extreme value of a (= 2m/3) or can be avioded by choosing multiple metric 
patches for nonextreme values of a, the junctions again having only C° conti- 
nuity. In either case, the tidal forces depend on controllable parameters m and 
[3 so that practical traversability is assured. As an illustrative seed solution, 
we considered only the Ellis I solution in the Einstein minimally coupled scalar 
field theory and finally mapped the extended solution into that of Brans-Dicke 
theory 

In connection with the static Brans-Dicke I solution, we recall a long stand- 
ing query by Visser and Hochberg [24] which has been the guiding motivation in 
Sees. II- V: "It would be interesting to know a little bit more about what this re- 
gion actually looks like, and to develop a better understanding of the physics on 
the other side [that is, across the naked singularity] of this class of Brans-Dicke 
wormholes. " The analyses above answer how one could achieve a both sided 
asymptotically flat traversable singularity free wormhole via Wick rotation of 
the Ellis I solution which is merely the conformally rescaled Brans-Dicke I solu- 
tion. Because of the fact that Ellis III is rcderivable from Ellis I and that the 
former is a traversable wormhole, we can say that Brans-Dicke I solution is also 
a traversable wormhole, but only in its Ellis III reincarnation. We believe that 
this argument provides some understanding of the "other side" : The mathemat- 
ical operations of conformal rescaling plus Wick rotation plus a trigonometric 
relation eases out the naked singularity and converts the Brans-Dicke I solution 
into a traversable wormhole. 

The next task was to investigate into the extended solutions with new pa- 
rameters. Thus, in section VI, we studied the geodesic motion in the extended 
geometry and obtained remarkable new results in the Einstein minimally cou- 
pled theory: For nonzero values of constant 9, the spacetime acts like a radial 
sink. For 9 = 0, the spacetime allows nonradial motions (U 3 ^ 0) but the 
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Matos-Nunez parameter a can be entirely expressed in terms of the constants of 
motion. In Sec. VII, we calculated the Sagnac effect in the extended spacetime 
and found that the delay depends on a. If a = 0, the delay is zero imply- 
ing that a could be interpreted as a rotation parameter thereby supporting the 
conclusion of Matos and Nunez [8] from an altogether different viewpoint. 

A very pertinent question is the stability of wormholes in the scalar-tensor 
theory, especially in the Einstein minimally coupled theory under consideration 
here. To our knowledge, the instability of black hole solutions in the Einstein 
conformally coupled theory was shown first by Bronnikov and Kireyev [23] . As 
to the wormholc solutions, it has been shown by Bronnikov and Grinyok [17] that 
they are also unstable under spherical perturbations Sip of the scalar field satis- 
fying tp : £ = 0. Instabilities occur due to the occurrence of negative energy levels 
in the effective potential. The authors suspected that such instabilities could 
be generic features of scalar-tensor theories. Instability of massive wormholes 
in the Einstein minimally coupled theory has been shown by Armcndariz-Pfcon 
[17]. However, he also showed that small mass wormholes in that theory could 
still be stable. It should be of interest to study whether the extended solutions 
are stable under perturbations. This remains a task for the future. 
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Appendix 

A wormhole is defined to be traversable by a hypothetical traveler if it satis- 
fies some general constraints [1]. We shall demonstrate that the Ellis III worm- 
hole satisfies all of these. It can be made traversable even by a human traveler 
under suitable choices of constants to and (3. Let us put one branch, say, the 
+ve branch of Eq.(29) into the metric (24) and rewrite it in the standard MTY 
form [1] by defining a radial variable p as 

{I 2 + to 2 ) exp[2/3{ J - arctan(-)}] = p 2 . (Al) 
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(Note that I — ► ±00 implies p — » ±00 and Z — » ±0 implies p ime*' 5 .) Then 
the metric (24) in the coordinates (t, p, 9, ip) becomes 



ds 2 = -e^'dt 2 + ^4(p) + P 2 ( d ° 2 + sin2 W 2 ) 



(A2) 



where the rcdshift function $ is 



$(p) - 



arctan 



\ m J 2 



and the shape function & is 

[J(p) - m/3] 2 



Hp) 



1 - 



exp[2/3{| -arctan(^)}] 



(A3) 



(A4) 



General constraints on b and <f> to produce a traversable wormhole are satined 
by the functions in (A3) and (A4). It may be verified that: (1) Throughout 
the spacetime, 1 — ^y- > and ^y- — > as p — > ±00 (2) Throat occurs at the 
minimum of p where &(p) = p. This minimum p corresponds to l — m(3 so 
that 

p = m(l + /3 2 )5 exp[/3{- - arctan/3}]. (A5) 

(3) The spacetime (A2) has no horizon, that is, $ is everywhere finite. (4) The 
coordinate time t measures proper time in asymptotically flat regions because 
$ — > as p — > ±00. (5) The spacetime has no singularities, as discussed in 
Sec.IV. 

Some more constraints are necessary if the trip is to be undertaken by a 
human traveler: (a) Trip begins and ends at stations located on cither side of 
the throat where the gravity field should be weak. This demands that (i) the 
geometry at stations must be nearly flat, or, <C 1, (ii) the gravitational 
rcdshift of signals sent from stations to infinity must be small, or, |$| <C 1 
and (iii) the acceleration of gravity at the stations must be less than one Earth 



gravity g§ — 980cm sec 2 , or, 



While the first two con- 



\ pi dp ~ 

straints (i) and (ii) are easily met in virtue of the general constraints (1) and 
(4) respectively, (iii) gives 



m(3 



l 2 +m 2 



< 



9® 



(A6) 



For fixed finite values of m and /?, e* — ► 1 for large I and hence this constraint 
can be easily satisfied at the stations, (b) The tidal forces suffered by the 
human traveler should be tolerable which means that the magnitude of the 
differential of four acceletation \A~ct\ should be less than g e in the orthonormal 
frame (eg, , e y , ej, eg, ) of the traveler. This constraint translates, for a traveler of 
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length (head to foot) <~ 2meters, into the following bounds on the components 
of curvature tensor computed in his/her frame: 



< 




dp 



(A7) 



x 2metr. (10 10 cm) 2 ' 



This bound is meant to constrain $(p) which is already well behaved at the 
throat p = p where its value is (3 [arctan/3 — §] . This value is —1 as f3 — ► oo, 
and as f3 — > 0. In general, as p — > ±oo, $ vanishes. Thus the constraint 
(A7) is easily satisfied. This result is expected since all the curvature tensor 
components fall off with distance and vanish at infinity [19]. Let us look at the 
lateral bounds given by 



iL 



< 



iL, 

I 3'0'3'0' 



9® 



2p 2 



db 
dp 



d$> 

2(p-b)^- 
dp 



(A8) 



c 2 x 2metr. (10 10 cm) 2 



which constrain the speed v with which the traveler crosses the wormhole, 7 = 
(1 — v 2 /c 2 )~2. For the metric (A2), the above works out to 



iL 



< 



iL, 

I 3'0'3'0' 



9® 



7 

2p 2 



(- 

V c 



2m(m + lf3) 
l 2 + m 2 



2m(3{l - m(3) 
I 2 + m 2 



J) 



c 2 x 2metr. (10 10 cm) 2 



The maximum tidal force should be experienced by the traveler at the throat 
at I = m/3 or equivalently at p = p Q . Hence the velocity v at the throat is 
determined by the following constraint 



Pi 



< 



1 



(10 10 cm) 2 ' 



(A10) 



It is possible to adjust m and [3 such that we have a throat radius p < 
(say). Assuming that v <C c or 7 ~ 1, we obtain v < 30metr./s(jQ^ 
shows that the speed v across the hole could be made reasonably small, (c) The 
traveler should feel less than g® acceleration throughout the trip which requires 



lOmetr. 
. ■) which 



1 - 



d_ 

dp 



( 7 e*) 



< 



(AH) 



With 7 ~ 1, this works out to the same constraint as in (A6) which is already 
satisfied, (d) The total proper time interval At measured by the traveler and 
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the coordinate time interval At measured at the stations for the entire trip 
should be about a year (say), that is, 



/ dL 
— < lyear (A12) 
- Ll «7 



/ 2 dT 
— < ^ear (A13) 
-Li 

where L = —L\ and L = +L2 arc the proper distances of the stations measured 
from the throat. Since 7 <~ 1, At ~ (L 2 + Li)/v where v is the average velocity 
not exceeding 30metr./s for a 10 metr. throat radius. This suggests that the 
total separation between the stations should be of the order of 9.4 x 10 5 Km. 

It is expected that the rotating wormhole should also be traversable, the 
only physical difference is that the travelers would feel an additional centrifugal 
force. However, the mathematical treatment of traversability criteria in the 
rotating case is quite cumbersome. We leave it as a future task. 

Consider the rotating solution, metric (7), with / generated from the +ve 
branch of Ellis III solution (29) given by 

8m/3(5exp[2/3{+f - arctan(-M}l 

f(l;m,(3,a) = 9 1 X 2 Vm , . 

M ' ' a 2 +4<5 2 exp[4/3{+f ~arctan(^)}] 

Since the metric components in (7) do not depend on t, the vector is a Killing 
vector of the spacetime but is not everywhere orthogonal to the constant t— 
hypcrsurfaces. Let us use a time variable T by 

dT = dt + a cos 9dip. 

The vector is not a Killing vector but is everywhere orthogonal to the con- 
stant T— hypersurfaces. Hence the metric (7) becomes 

ds 2 = -fdT 2 + r 1 [dl 2 + (I 2 + l 2 ){d0 2 + sin 2 9d^ 2 )] . 
Under the redefinition dT = dt' + dl, the metric (A) becomes 



ds 2 = -dt' 2 + (dl - v/1 - fdt'f + (I 2 + ll){d0 2 + sin 2 0dip 2 ). 

It was shown in Eq.(63) of Sec. VI that one solution of the geodesic equations 
was radial motion (U 3 = 0) on the plane 9 — const. Without loss of generality 
we can choose 9 = ir/2 which implies that dT — dt on the equatorial plane. The 
metric on the constant t' surface for 9 = n/2 is then da 2 = dl 2 + (I 2 + l 2 )dip 2 . 
This surface is isometric to a catenoid {(x, y, z) \ (x 2 + y 2 )? — Z cosh(p)} in 
E 3 , the radius of the central hole being lo [19]. The catenoid with a central hole 
is the sink referred to in Sec. VI. The generators of the catenoid are the radial 
trajectories of a hypothetical traveler. 
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